Instanton Approach to Josephson Tunneling between Trapped Condensates 



Yunbo Zhang 1 ' 2 * and H.J.W. Miiller-Kirsten 1 * 

1 Department of Physics, University of Kaiserslautern, D-67653 Kaiserslautern, Germany 
2 Department of Physics and Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, P.R. China 

" l 1 An instanton method is proposed to investigate the quantum tunneling between two weakly- 

^ , linked Bose-Einstein condensates confined in double-well potential traps. We point out some in- 

. trinsic pathologies in the earlier treatments of other authors and make an effort to go beyond these 

' very simple zero order models. The tunneling amplitude may be calculated in the Thomas-Fermi 
approximation and beyond it; we find it depends on the number of the trapped atoms, through the 

. chemical potential. Some suggestions are given for the observation of the Josephson oscillation and 

/"N ' the MQST. 

^ ' PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj 
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"I I. INTRODUCTION 

13 ; 

The first experimental observations of Bose-Einstein condensation (BEC) in dilute gases of trapped alkali atomsQ| 
have stimulated the studies of condensates in double well traps. With a fascinating possibility of the observation 
of new kinds of macroscopic quantum phenomena which are related to the superfluid nature of the condensates, 
H ■ numerous authors have addressed this subject area both experimentally and theoretically. Some recent experiments 
, t \ have investigated the relative phase of two overlapping condensates in different hyperfine states]^] and robust inter- 
ference fringes between two freely expanding condensates have been observed j| after switching off the double-well 
potential that confines them, indicating phase coherence both in space and time. In fact the possibility of condensate 
tunneling between two adjacent atomic traps and detection of Josephson-like current phase effects have been previ- 
ously suggested [|| 0, |?| and intensively studied by two main approaches which are capable of dealing with quantum 
tunneling in this variant of the Josephson effect. Authors with a quantum o ptic s background tend to favor models in 
q ■ which two boson modes are involved, the so-called two-mode model||, [| |l(], jy], [l^, [l3|, |l4], The other category of 
theories is based on using the differences of condensate phases and atom numbers between the two sides of the trap 
as conjugate quantum variables^]. A detailed comparison between these two approaches is given in ref . []r7f . 

The superfluid nature of condensates can be fully tested only through the observation of superflows. Despite the 
^-j- . many experimental efforts being focused on the creation of a Josephson junction between two condensate bulks, 
' direct experimental evidences for the atomic oscillation are still far from being realized. Theoretically the Josephson 
junction problem has been studied in the limit of noninteracting atoms for small- amplitude Josephson oscillations (|, ||, 
■ including finite-temperature (damping) effec ts |lq| . Dynamic processes of splitting a condensate by raising a potential 
barrier in the center of a harmonic trap |l7], |l8|| and decoherence effects and quantum corrections to the semiclassical 
mean- field dynamics || O, [ll], have also been studied. It has been pointed out|| that even though the 

Bose Josephson Junction(BJJ) is a neutral-atom system, it can still display the typical ac and dc Josephson effects 
occurring in charged Cooper-pair superconducting junctions. Moreover, a novel nonlinear effect has been predicted 
to occur in this BJJ: The self-trapping of a BEC population imbalance arises because of the interatomic nonlinear 
interaction in the Bose gasj|, This was considered to be a novel "macroscopic quantum self-trapping" (MQST) 
and was predicted to be observable under certain experimental conditions. Three related parameters, i.e. the ground 
state energy E°, the interaction energy U, and more importantly, the tunneling amplitude K, are still undetermined 
for a specific geometry of the potential well and have been taken as constants in refs.j8|, |9|. It is the main purpose of 
this paper to present a rigorous derivation of these quantities and we find that they actually depend on the number 
of atoms TV. This A^-dependence refines the conclusions and makes the self-trapping easier to observe. 

We develop the instanton method for a sensitive and precise investigation of the tunneling between two condensates. 
The almost trivially looking problem of the tunneling behavior in a double-well potential has attracted much attention 
from theorists for decades. For a single particle, the solution can be found even in quantum mechanics textbooks fp2|. 
The advantage of a nonperturbative method, as presented here, is that it gives not only a more accurate description of 
the tunneling phenomena but also a comprehensive physical understanding in the context of quantum field theory. The 
periodic instanton configurations, which have been shown to be a useful tool in several areas of research such as spin 
tunneling p3[ p4j, bubble nuclcationpa] and string theory pq], enable also the investigation of the finite temperature 
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behavior of these systems. In the case of the Bose-Einstein system, however, we need to evaluate the tunneling 
frequency for a finite chemical potential even at zero temperature, due to the nonlinear interaction between the 
confined atoms. Therefore the chemical potential here replaces the position of the excited energy and gives rise to an 
expected higher tunneling frequency. 

Thus, in this paper we will study the atomic tunneling between two condensates in a double-well trap, with special 
attentions paid to the N— dependent tunneling amplitude. In Sec. II we review the two mode approach and summarize 
3 critical conditions which are related to the MQST effect. Analytical solutions to the coupled nonlinear BJJ equations 
are derived in Sec. Ill and some novel features about the different modes of MQST are discussed in detail. In Sec. 
IV the often used Thomas-Fermi approximation and the corrections beyond it are used to obtain the two important 
parameters in the BJJ tunneling system. By means of the periodic instanton method, wc calculate the tunneling 
amplitude/frequency between the two condensates in Sec. V and find they actually depend on the number of the 
trapped atoms. Our results can be compared to the simple single-particle tunneling result in ref.Q and we find that 
the latter corresponds to the low-energy or noninteracting case. A detailed discussion about the optimal condition 
for observation of the MQST and the Josephson tunneling is the subject of Sec. VI where some comments are made 
to the maximum amplitude of the oscillation. Finally our summary and conclusions are given in Sec. VII. 



II. TWO-MODE MODEL AND MACROSCOPIC QUANTUM SELF-TRAPPING 

It was realized long ago that the ground state properties of a trapped Bose-gas can be well described by the 
Gross-Pitaevskii equation(GPE) |27| 



9$ k 2 , 
dt 2m 



V ext (r) + g \<S>\ 2 $ (1) 



where V ex t is the external trapping potential and go — i7Th - is the interatomic coupling constant, with a, m the 
atomic scattering length and mass, respectively. In this paper we will consider systems interacting with repulsive 
forces with a > 0. To obtain the ground state properties, one can write the macroscopic wave function(or from the 
viewpoint of phase transitions, the order parameter) of the condensate as $(r, t) — <j)(r)e~ ltJ ' t / h , where a condensate 
in the stationary state was assumed. Then the GPE (Q) becomes 

fi 2 V 2 

Mr) = (H + . 9o 2 (r)) 0(r), H = -— - + V ext (r) (2) 

2m 

with Ho the Hamiltonian for the condensates of noninteracting bosons. This equation explains the significance of 
the chemical potential \x as an energy of the stationary level and can be used further for calculation of the coupling 
between two condensates, say, the tunneling amplitude(transfer matrix element or Josephson tunneling term) as will 
be shown below. 

Consider a double-well trap produced, for example, by a far-off-resonance laser barrier that cuts a single trapped 
condensate into two parts. In the barrier region the modulus of the order parameter in the GPE is exponentially 
small. In other words, the overlap between the condensates occurs only in the classically forbidden region, where the 
wave function is small and nonlinear effects due to interactions can be ignored. Thus we look for the solution of the 
time-dependent GPE (j|) with the two-mode variational ansatzg, 

$(r, t) = Mt)^i(r) + lfe(i)*2(r). (3) 



Here, the complex coefficients ij)i(t) = \J iVj (t) exp [i0i (t)] are spatially uniform and contain all of the time dependence, 
while the two states $i(r) and $2(1") are localized in the left and right wells, respectively, and contain all of the 
position dependence. The total number of atoms is conserved in the macroscopic quantum tunneling process, i.e., 
N% + N 2 = \ipi\ 2 + iV'al = Nt- In a more accurate theory, one could account for the slow time evolution of the states 
$i(r) due to the mean field interaction since both the number of particles Ni(t) and phases in each well evolve 
in time. This, however, can be shown to have a negligible effect when the amplitude of oscillation is relatively small 

©• _ _ 

Substituting the two- mode ansatz (||) into the GPE (|l|), multiplying by $i(r) and integrating over position wc 

obtain the equations of motion for the complex coefficients ^>i(t), i.e. the BJJ equations 

ih^ = (Ex + UxNx^i-Kih, (4) 
ih^ = {E 2 + U 2 N 2 )^ 2 -K^, 
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where damping and finite temperature effects are ignored. Here £1,2 are the zero-point energies in each well, 

E h2 = J dr$ 1)2 (r)ffo$i >2 (r), (5) 
UiNi are proportional to the atomic self-interaction energies, with 

UL2 = .90 J dc*\ t2 (T) (6) 
and K describes the amplitude of the tunneling between condensates 

K = - [ Wi(r)F„$ 2 (r). (7) 



These quantities, which are expressed in terms of appropriate overlap integrals of the wave-functions $1.2(1"), have 
been taken as constants in the previous works and it is one of the main tasks of this paper to obtain them analytically, 
provided that a specific geometry of the trap is given. We further note that the (real) ground state solutions $1,2(1") 
for isolated traps with equal population in each well, i.e. N\ — N2 — At/2, satisfy the orthonormality condition 



Jdr$ i (r)$ j (r) = 5 ij . (8) 



In the derivation above we neglected terms such as J Q\&2<lr, / $i$2rfr, and J $i$2^ r which have the meanings 
of higher order overlaps between the condensates. Attempts to include the mean-field contribution to the coupling 
between the condensates have also been made|2^]. In refs.^O) the authors studied the coherent atomic oscillations and 
resonances between coupled BECs with time-dependent symmetric and asymmetric trapping potential and oscillating 
atomic scattering length and an interesting chaotic macroscopic quantum tunneling phenomenon was found to exist. 
Our analytical result for E, U and K may be useful in the observation of this macroscopic quantum effect. 

Straightforwardly we can obtain the equations for the relative phase <j>(f) = 62(f) ~ 61(f) and fractional population 
imbalance z = [Ni(t) - N 2 (t)] /N T §, [fif 



z = — y/l — z 2 sin 4>, 

6 = Az H — -=^= cos 4> + AjE, (9) 
V 1 — z 2 

where we have rescaled 2Kt/h to a dimensionless time t for time- independent coupling K between two condensates. 
The parameters AE and A, which determines the dynamic regimes of the BEC atomic tunneling, can be expressed as 

= E 1 -E 1 Ui^U2 N 
2K 4K r ' 

A = U=(U 1 + U 2 )/2. (10) 

The canonically conjugated variables z and 0, satisfying the corresponding Hamilton canonical equations of motion, 
z = — ^ and <j> = 0^, suggest that the total energy of the above system is conserved, 

A 



H = —z 2 - VI - z 2 cos <f> + AEz. 



In a simple mechanical analogy, H describes a nonrigid pendulum of tilt angle <p and length proportional to \/T — z 2 , 
that decreases with the angular momentum z. For simplicity we restrict our calculations to a symmetric double- well 
potential, which means that the spatially averaged quantities (^) and (|^) are equal for each well, allowing us to make 
the simplifications AE = 0. In this case, the equations of motion reduce to 



with the conserved energy 



z(t) = --yji - z(t) 2 sin 4>(t), 
z(t) 

<t>(t) - A z ( t ) + ^=A±= C0S cj>(t), (11) 

v 1 - z \fr 



H = H[z(t), cf>(t)] = ^z(tf - s/l-z(t) 2 cos 4>(t) 

ff[z(O),0(O)] = ^z(Q) 2 - ^1-^(0)^08 0(0) = H a . (12) 
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In order to see the oscillations of the fractional population imbalance and the different kinds of (running- and 
7r-phase) MQST modes more transparently, we introduce an effective classical particle whose coordinate is z, moving 
in a potential W(z) with the initial energy Wq. The effective equation of motion is 

z{tf + W(z) = Wo (13) 



where 



W(z) = z 2 



A 2 



(14) 



W = 1 - Hi = W[z(0)} + i(0) 2 . (15) 

Increasing the value of 1 — AHo from negative to positive changes the effective potential W(z) from a double-well 
to a parabolic. The motion in the parabolic potential is Rabi-like oscillation with a zero time-average value of the 
fractional population imbalance z. For fixed parameters A and Hq, the oscillations with small effective energies 
Wo are sinusoidal, the increasing of the effective energies adds higher harmonics to the sinusoidal oscillations. In 

the case of double- well potential with two minima located at z± = ± v A — - and the barrier height given by 

W(0) - W(z±) = (Ag °7 1)2 , the motion is very different from that in the parabolic potential. If the effective energies 
is larger than the barrier between two wells, that is, Wo > or Hq < 1, the motion is a nonlinear Rabi oscillation with 
a zero time-average value of z, which corresponds to the periodic flux of atoms from one condensate to the other. If 
the effective energy is lower than the potential barrier, Wo < or Ho > 1, the particle is forced to become localized 
in one of the two wells and the population in each trap oscillates around a nonzero averaged (z(t)) ^ which has 
been termed MQST in |], [l(J. Based on these analyses, we know the critical condition for MQST corresponds to the 
point where the effective energy equals the potential barrier, Wo = or Ho = 1. 

The lowest four stationary state solutions (i = 4> = 0) of eq.(pT|) are: a) the symmetric ground state (j> s = 2rm, 
z s = with energy E + = — 1; b) the antisymmetric eigenfunction <p s — (2n + 1)tt, z s — with energy _E_ = 1; and 

c) the z— symmetry breaking states cj> a — (2n + V)ix, z s = ±yl— with energy — | (A + j^j > 1, respectively. 
The degenerate eigenstates that break the z symmetry c) are obviously the result of the nonlinear interatomic 
interaction. In this sense we may verify again the MQST condition Ho > 1, which means the populations become 
macroscopically self-trapped with (z) ^ 0. From above the 2 lowest-lying states a) and b) with z 8 = are obviously 
not MQST(_ffo = —1,1), only the symmetry breaking states c) with z s ^ can be regarded as being self-trapped 

( ff o = H A + x)> 1 )- 
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FIG. 1: The parameter range for the appearance of MQST for 0— phase mode: MQST, J— > Josephson Tunneling. 
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In summary we observe 3 different critical conditions in refs. |9L 101] . One is the MQST condition which is defined 



< z{t) >^ => H Q > 1 
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|z 2 (0)- Vl-^(0)cos[^(0)]>l. 



(16) 



In a series of experiments in which 0(0) and z(0) are kept constant but A is varied by changing the geometry or the 
total number of condensate atoms, we have the critical parameter 



A, 



1 + yjl - 2(0) 2 COS 0(0) 

z(0) 2 /2 



(17) 



and the critical condition for MQST is A > A c . On the other hand, changing z(0) and keeping (f)(0) and A constants 
we have a critical initial population imbalance z c , which takes the value 



A 



VA - 1 



(18) 



in both 0— phase mode and 7r— phase mode case, which describe the tunneling dynamics with the initial or time- 
averaged value of the phase across the junction being and 7r, respectively. But for 0(0) = 0, if z(0) > z c , MQST 
sets in, and for 0(0) = 7r, z(0) < z c marks the region of MQST. In Fig. 1 we show the parameter range for the 
appearance of MQST and Joscphson tunneling for 0— phase mode, while for the 7r— phase mode this phase diagram 
manifests much richer structure as shown in Fig. 2. It should be pointed out that the vertical axis A can extend to 
much larger values, but we illustrate in the figure only the area of interest. For A > 2 and 0(0) = n, the effective 
potential W(z) always has a double-well structure and the system is self trapped for all values of z(0). We conclude 
that in the case of the 7r— phase mode the MQST will set in for most of the areas, leaving only a small area for actual 
Josephson tunneling. 

Fig. 2 
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FIG. 2: Different MQST mode for 7r-phase mode: J -* Josephson Tunneling, I -> MQST Type I, II -» MQST Type II, R -» 
MQST unbounded running phase mode. The vertical and horizontal lines correspond to the case of Fig. 7 and Fig. 5 of ref. 
[ p"o"| |, respectively. 

The second condition comes from the stationary z— symmetry-breaking solution c); the value of z s is just the 
boundary of the first- and the second-type of trapped states which differ by the time-average value of the population 
imbalance (z) . We enter into some detail of the ir phase mode. According to jio) , there are two kinds of such 7r-phase 
modes with MQST: 



7r-phase modes MQST Typel <^ (z) < z s ^ 
vr-phase modes MQST Typell <=> (z) > z s ^ 

with z s being the stationary z-symmetry breaking value of the BJJ equations (^). Once A exceeds the value 

1 



A., 



(19) 



(20) 
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(cf. z s ~ ±y 1 — -^j-), a changeover occurs at the stationary state and the system goes from Type I of the 7r-phase 

trapped state into Type II. In the phase plane portrait of the dynamical variables z and (Fig- 7 in ref . |l0| ) for 
the fixed value of z(0) = 0.6, this critical condition corresponds to the 4 point-like trajectories located at (±tt, ±0.6). 
However, in consideration of the case of a fixed value of A, for example, Figs. 5c and 5d in ref. |u|, we find that this 
definition of the Type I and Type II 7r— phase modes will lead to a dilemma in the analysis of the phase diagrams and 
modification to this definition becomes necessary. We return to this point in the following section. 

The third condition defines whether the phase <p(t) is an unbounded running mode or is localized around tt. 
Again from the z ~ <j> phase diagram we know this boundary line can be determined as 



TT d(j) A 

2' ^ U=±1 =°^2 



U=±i = -, — l»=±i = o = #0 (21) 



which gives the critical values 



A3 = 7T^Wy or ^V 1 "^ (22 > 

In Fig. 2, we managed to give a phase diagram for the tt— phase mode, with special attention paid to different MQST 
modes. Here the area C shows again the Josephson tunneling area, others are 3 MQST areas, with area I the bounded 
MQST of Type I, area II the Type II, and Running Mode for the unbounded running phase mode, respectively. We 
used the formulae (|l7|), ( p0[ ) and (^2|) for these three curves, from the bottom up, respectively. We point out that the 
second and the third conditions apply only to the (j>(0) = tt case. 



III. EXACT ANALYTIC SOLUTIONS AND DIFFERENT MQST MODES 



From the equation of motion for the effective classical particle ([L3|) we can easily obtain the real time exact solutions 
in terms of the Jacobian elliptic functions cn and dn. This has been done in refs. [ll| but there are some notation 
errors in the corresponding expressions. We thus rewrite the solutions here and present a detailed discussion, especially 
on some featured oscillation modes. For the parameters in the area of Josephson tunneling the population imbalance 
oscillates sinusoidally or nonsinusoidally between a positive initial value C and a negative one — C, 

z(t)=Ccn(^(t-t ),kj (23) 



where 



c2 = 3s((*oA-i) + £ 

c? = ±(^-(H A-l)) (24, 
C 2 = 2 VA 2 + 1 - 2H A 



with < k < 1 the elliptic modulus, 



2 

C 2 2 V C ) 2 ^ ' VA 2 + 1 - 2# A 



In the case of MQST, however, the elliptic function cn will be replaced by its counterpart dn and the oscillation is 
about a nonzero average value (z) 

z(t)=Cdn(^(t-t ),l/kY (26) 

From this we know the value z(t) — is inaccessible at any time and the modulus is now 1/fc with k > 1. The 
integration constant £o can be determined from z(t)\ t =o = z(0) : 

2k , fz(0) ,\ 2 / 2(0) ,\ 

*o = 77*- cn " )=^^^ F arccos-^/c (27) 



CA \ C 1 J AV« 2 + C 2 V C 
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where F ((f), k) = Jft d(f>(l — fc 2 sin</>) 1 I 2 is the incomplete elliptic integral of the first kind. Noting the following 
correspondence 

< k < 1, a 2 > 0, Ho < 1, Josephson Oscillation (28) 
k > I, a 2 < 0,H a > 1, MQST 

we easily derive the physical condition for the onset of MQST, i.e., Ho — 1 and A = A c , from the mathematical 
condition k = 1, where the character of the elliptic function solution changes. The Jacobian elliptic functions cn(u, k) 
and dn(u, k) are periodic in the argument u with period 4JC(k) and 2/C(fc), respectively, where K,[k) is the complete 
elliptic integral of the first kind. The time period of the oscillation of z(t) is then given by 

oh 

t = _4£(fc), for < k < 1, (29) 
t = -|-2/C(l/fc), for k > 1. (30) 



In the limit of small amplitude, or linear oscillations, k — ► 0, K(k) — ► tt/2, we have 

2tt 

^/A 2 + l-2i?oA = ' 

For 0,7r— phase oscillations, Hq — > =pl, the periods in scaled units 2Kt/T% become 

2n 



vO±A 

respectively, which agree with the unsealed ones (eqs (4.7) and (4.10) in ref . [|To| ) 



(31) 



(32) 



r -* = ^/4K 2 ± 2UN T K/2Trh. (33) 

We consider here two special cases, i.e., the 0,7r— phase modes. For the 0— phase mode, by inserting cj>(0) = into 
eq.(EJ), we may show that the oscillation amplitude of z(t) is just the initial value z(0), 



C = z(0). (34) 

On the other hand, for the 7r— phase mode, this amplitude will depend on the values of A and z(0), 

C 2 = z 2 (0), when A < A s or z(0) > z„, (35) 

4 

A 2 



C 2 = z 2 {0) + — AVi - z 2 {0) - 1 , when A > A s or z(0) < z s 



In the case of Josephson tunneling for both and 7r— phase modes, that is, when the oscillation is in the form of 
elliptic function cn-type, we always have C = z(§). The reason is that under the condition A > A s or z(0) < z s for 
it— phase mode the system is obviously in the region of MQST type II. Therefore we should use the elliptic function dn 
for MQST instead. Generally for MQST the average population imbalance may be calculated as(the elliptic function 
dn oscillates between 1 and k' dn ) 



<-> = ^(l + ^) = ^(l + Vl _ W)-f l + V 1 -^ V ( 36 ) 



For the 0— phase mode, C — z(0) and ( 2 = 2 ( A-\/l — z 2 (0) + 1), the average population imbalance is 

z(0) 



4(A v /l-z 2 (0) + l) 
1 + \ 1 z(0) 2 A 2 



(37) 



We take the values of Fig. 2 of ref.JT(| as an example, A — 10, z(0) ~ 0.65, (z) — 0.46511, which is exactly the case. 
For the it— phase mode, there are two cases: 
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Case I: A < A s or z(0) > z s we have C — z(0) and Q 2 — 2 — Ky/1 — z 2 (0)Jwhich give the average population 
imbalance as 



(z) 



z(0) 



1 + \ 



4 (l - AVl-^(O)) 
1 z(0) 2 A 2 



(38) 



For the parameter in Fig. 3d of ref.jlOt, z(0) = 0.6, A = 1.2, (z) = 0.54944 < 0.6. This is again exactly the case. 
Case II: A > A s or z(0) < z s we have 



4 r 



C 2 = z 2 (0) + — Ay/1 - z 2 (0) - 1 
A z L 

C 2 - 2(a v /i-^ 2 (o)-i), 



(39) 



the average population imbalance is the same as that in Case I eq. (|3q) . For Fig. 3f of ref.jllj, z(0) = 0.6, A = 
1.3, (z) = 0.63715 > 0.6. 

We may compare our phase diagram Fig. 2 with those previously known results. First we have a look at the Fig. 
7 of ref . ]l0| , which illustrated the behavior of the system in the z — 4> phase-plane, and concentrating only on the 
tt— phase mode. At a fixed value of z(0) = 0.6, increasing the value of A will bring the system from area C to area I 
(Type I bounded MQST with (z) < 0.6), then area II(Type II bounded MQST with (z) > 0.6), and finally into the 
running mode area. This process is described by a vertical line in Fig. 2 and it intersects successively the three curves 
at 



A c = 10/9 = l.i, A, = 1.25, A 3 = 2.5, 



(40) 



from the bottom up. In Table I we list detailed descriptions and graphical features of the oscillations in accordance 
with Fig. 7 of ref. pQ| . The coincidence here is obvious. 



Table I: Different oscillation modes for fixed value z(0) = 0.6 



(f)(0) = TT 


Descriptions 


Graphics 


A = 


Rabi oscillation around (z) = 


O 


< A < A c 


sinusoidal/nonsinusoidal oscillation around (z) = 





A = A c = 1.111 


the trajectory shrinks and is pinched at the point z = 
marking the onset of 7r-phase MQST type I with (z) < z(0) 
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A c < A < A s 


upon further increasing the area are divided into 2 parts 


O 
O 


A = A s = 1.25 


the 2 divided areas collapse to 2 point-like trajectories 
at z = 0.6; boundary between MQST type I and type II 


1 


A s < A < A 3 


further increase of A induces a reflection of the trajectory 
about the fixed point and 7r-phase MQST type II with (z) > z(0) 


o 
o 


A = A 3 = 2.5 


boundary between bounded-mode(type II) and running-mode MQST 


V 

A 


A > A 3 


the trajectory joins the unbounded running-mode MQST 


UJ 
777. 



Secondly we go through Fig. 5 of ref. |X(J to sec what happens if we keep the value of A constant and increase 2(0). 
This is indicated in Fig. 2 by a horizontal line. For small values of z(0), the phase is unbounded and the system 
exhibits running phase MQST. However, above a certain value z 3 = y/T— 4/A 2 (cf. A 3 ; there is an error in ref. ficfl 
where the condition z(0) = 2z s = 2y/l — 1/A 2 is unphysical since it is larger than 1), with (z) still nonzero, the 
phase becomes localized around tt and remains bounded for all larger values of z(0). Then the puzzle appears here. 
According to ref. 0, in their Figs. 5(c) and 5(d), z(0) = 0.7 and 0.98 mark the two different types of tt— phase MQST 
since they are on either side of the stationary state value of z s = y/l — 1/A 2 . In this case, along the direction of the 
horizontal line from left to right, we should enter Type I first (z(0) = 0.7 < 0.92), then Type II(z(0) = 0.98 > 0.92). 
But this is not the case, instead we enter Type II first. The horizontal line cuts the two curves at z 3 = 0.6 and 
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z s = 0.92, with the latter corresponding to a point-like trajectory in Fig. 5d of ref. |M at z = 0.92, <f> = it. Another 
possibility is that although the initial value z(0) = 0.7 < 0.92 but the average value (z) > 0.92. However, from eq. 
(|38|), we find that in both cases, i.e. for z(0) = 0.7 and 0.98, we have (z) < 0.92. 
A solution to this dilemma is to define the MQST types I and 11 as follows: 

MQST Type I & (z) < z(0) (41) 
MQST Type II ^ (z) > z(0) 



instead of (19). Under this novel definition we can explain the phase diagrams as follows: For z(0) = 0.7 the average 
population difference (z) > 0.7 which is obviously Type II, while for z(0) = 0.98 we have (z) < 0.98 which is type 

I according to this definition. Therefore, along the direction of the horizontal line from left to right, we enter Type 

II first ((z) > z(0) = 0.7), then Type I((z) < z(0) — 0.98). At the threshold point, the average value coincides with 
the initial value, (z) = z(0) = z s = 0.92, and the trajectory shrinks to a point. We also notice that the oscillation 
trajectories for z(0) = 0.8 and 0.98 are the same. However, they belong to different regions of MQST. In Table II 
we illustrate the oscillation behavior of a 7r— phase mode when the value of z(0) is increased while keeping A = 2.5, 
which can be re gard ed as an additional remark to Figs. 5c and 5d in ref. JTo|] . A further evidence of the correctness 
of the criterion ( |4l|) comes from the fact that from the analytical result for the average population imbalance in the 
it— phase mode, that is, eq. d38|), we can derive the second critical condition for A, 

(*) - z(0) =► A, = 1 (42) 

y/1 - Z 2 (0) 

which agrees with eq. (pfj|) . And in the case of the 0— phase mode there is no solution for (z) = z(0) with (z) determined 
by eq.(|37|), indicating that the self-trapping can only appear as MQST type I((z) < z(0)). 



Table II: The oscillation behavior of a ir phase mode for fixed value A = 2.5 



z(0) 


The Potential W(z) 


The Trajectory </> 


Descriptions 





double-well with Wo 
just on the top of the barrier 


pinched 
at z(0) = 


always 
MQST 


< z(0) < z 3 


Wo under the barrier whose 
height increases with z(0) 


unbounded 


running-phase 
MQST 


z(0) = z 3 = 0.6 


barrier height increases 
continuously with z(0) 


0U=l=7T/2, 
dz l*=l - U 


boundary between 
unbounded and localized <p 


Z3 < z(0) < z s 


barrier height increases 
continuously with z(0) 


localized 
around n 


bounded MQST 
type II with (z) > z(Q) 


z(0) = z s = 0.92 


barrier attains 
its maximum height 


shrinks to point-like 
trajectory at 7r 


boundary between 
MQST type I and type II 


z s < z(0) < 1 


barrier height decreases 
with increasing z(0) 


expands but 
again localized 


bounded MQST 
type I with (z) < z(0) 


z(0) = 1 


barrier height decreases to the 
same height as z(0) = 0.6 


agrees 
with z(0) = 0.6 


bounded MQST 
type I with (z) < z(0) 



IV. THOMAS-FERMI APPROXIMATION AND BEYOND 



The macroscopic wave function $ associated with the ground state of a dilute Bose gas confined in the potential 
Vext (r) obeys the GPE (|l|) , which can be obtained alternatively using a variational procedure: 
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where the energy functional E is defined by 

= J d 3 r 

= Ekin+Eho + Eint. (44) 



^|V$| 2 + Wr)|$| 2 + f |<I>| 4 
2m 2 



The three terms in the integral are the kinetic energy of the condensate Ekin, the (an)harmonic potential energy Eho, 
and the mean-field interaction energy Ei nt , respectively. By direct integration of the GPE (||) one finds the useful 
expression 

Nfi = E km + E ho + 2E mt (45) 

for the chemical potential in terms of the different contributions to the energy functional (|44|). Further important 
relationships can be found by means of the virial theorem 

2-Efcm — 2Eho + 3Ei n t = (46) 

and the thermodynamic definition of the chemical potential 

fj, = dE/dN. (47) 

It was shown |^|, [3l[] that the ground state properties of the trapped Bose gas with repulsive interactions may be 
described quite accurately for a sufficiently large number of particles N by a Thomas-Fermi approximation(TFA), in 
which the kinetic energy, which is also usually named quantum pressure, is neglected completely. One then gets the 
density profile in the form 

n(r) = ^ 2 (r)= 5o - 1 [M-Kxt(r)] (48) 

in the region where \i > V ex t( r ), and n(r) = outside. In the simplest case of an isotropic harmonic trap V ex t{r) = 
towqT 2 /2, the normalization condition on n(r) provides the relation between chemical potential and number of particles 



flTF 



huJo ( 15Na\ 



2/5 



(49) 

V aho J 



where the harmonic oscillator length a; lo = (h/mtuo) 1 / 2 is introduced for simplicity. The density profile (E§) has the 
form of an inverted parabola, which vanishes at the classical turning point R defined by the condition fiTF — V ex t(R). 
For a spherical trap, this implies htf — mui 2 R 2 /2 and using result (|49| ) for /itf, one finds the following expression 
for the radius of the condensate 

/l57Va\ 1/5 

R = a ho (50) 

V ato ) 

which grows with N. The energy components assume the following simple values in the TFA 

Ekin n Eho 3 Ei n t 2 

^\r =0 ' — = r TF > — = r TF - (51) 

The above results are obtained by neglecting the kinetic energy term Ekin in the GPE and provide an accurate 
description of the exact solution in the interior of the atomic cloud where the gradients of the wave function are small. 
This approximation does not, however, take properly into account the decay of the wave function near the outer edge 
of the cloud, and the Thomas- Fermi wave function ( |4^ ) consequently leads to unphysical behavior for some properties, 
most notably the logarithmically divergence of the kinetic energy arising at the turning point. A good approximation 
for the density in the region close to the classical turning point can be obtained by a suitable expansion of GPE (||). 
In fact, when \r — R\ -C R, the trapping potential V ex t{r) can be replaced by a linear ramp, mu;gi?(r — R), and the 
GPE takes a universal formal yielding the rounding of the surface profile. Using this procedure it is possible to 
calculate the kinetic energy which, in the case of a spherical trap, is found to follow the asymptotic law[[| 

E "" V. r 4t4t^V (52) 



N 2 ' mR 2 \1.3a ho 
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Analogous expansions can be derived for the harmonic potential energy Eho and interaction energy Ei nt in the same 
large N limit |32|| . A straightforward derivation is obtained by using nontrivial relationships among the various energy 
components. A first relation is given by the virial theorem (fl6]). A second one is obtained by using expressions J45| ) 
and ([It]) for chemical potential. These two relationships, together with the asymptotic law (^2|) for the kinetic energy, 
allow one to obtain the expansions 

Eho 3 Ei n t 2 

= -jVTF + C, = —fJ-TF - C. (53) 

Correspondingly the chemical potential in TFA is modified beyond TFA as 

A* = A*tf + \c (54) 

We may now first calculate the two relevant quantities in the Bose Josephson Junction, i.e., the zero-point energy in 
each well E\ 2 in eq.(||) and the atomic self-interaction energies U\, 2 in eq.(^) both in TFA and beyond it, leaving the 
third one, the tunneling amplitude K , to the next section. We note that in the derivation above, the wave function 
is normalized to N . If one uses instead a wave function normalized to unity, the following correspondence should be 
realized 

U ia N ia - 2^, E ia ^ + (55) 

Setting the initial population difference z(0) = as in ref. |jj,i.e., Ni = N 2 = N, we obtain the ground state energy 
£1,2 and the interaction self energy Ui t2 Ni j2 for the isolated traps beyond the TFA 

E1.2 = ^tf + \c, UxflNip = ^tf - 2C, (56) 

while for TFA the correction terms disappear and these energies take simpler forms 

3 4 
Ei,2 = jfJ-TF, Ui j2 N 1}2 = ~htf- (57) 

Considering a condensate of N = 5000 sodium atoms confined in a symmetric spherical trap with frequency ujq = 
100Hz, we have the results beyond the TFA, E 1>2 = 1.18nK, U h2 N ia = 1.03nK, and in TFA E h2 = 0.9nK, Ui, 2 N lt2 = 
1.19nK, quite close to the values estimated in ref. [9. 



V. CALCULATION OF THE TUNNELING AMPLITUDE BETWEEN TWO TRAPPED 

CONDENSATES 



In this section we give a rigorous derivation of the amplitude of the atomic tunneling at zero temperature between 
two nonideal, weakly linked condensates in a double well trap. This induces a coherent, oscillating flux of atoms 
between wells, that is a signature of the macroscopic superposition of states in which the condensates evolve. To date 
there have been no reports of experimental observations of Josephson tunneling of a condensate in a double well trap. 
Josephson tunneling of a condensate in a one dimensional optical lattice was reported in ref. f3~if| and by controlling 
relative strengths of the tunneling rate between traps and atom-atom interactions within each trap, this technique 
has been used to produce atom-number-squeezed states in this lattice potential j35|. The analogy of the tunneling 
mechanism in a two-well potential and an array of wells makes it important to calculate the tunneling amplitude 
explicitly. 

The quantity K in eq. (|7j) corresponds to the tunneling amplitude which can be calculated by different methods, and 
we demonstrate in this work the use of the nonperturbative instanton approach. The nonlinear interaction between 
the atoms in the same well will be included, which modifies only the chemical potential fi both in and beyond the 
TFA. It is easily shown that this tunneling amplitude is just the quantity 1Z of B (up to a minus sign), if one observes 
the orthogonality property of the eigenfunctions eq.(||), with &i y2 (r) the local modes in each well. We study the 
amplitude for tunneling between the two condensates confined in the wells of an external double-well potential in the 
3-dimensional 

V ext {v) = ^rOc 2 - xl? + \^ly 2 + \muolz 2 (58) 
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where we have assumed as in ref . |g] that the interwell coupling occurs only along x and the wave function components 
in the other two dimensions are orthogonal and contribute a factor of unity. The two minima are located at ±Xq on 
the x— axis and the parabolic approximation to the potential in the vicinity of each minimum is 

V {2) (r - r 1>2 ) = X -mul [(x ± x ) 2 + y 2 + z 2 ] . (59) 
The barrier height between the two wells 

V = -mulxl (60) 

is assumed to be high enough so that the overlap between the wave functions relative to the two traps occurs only 
in the classically forbidden region where interaction can be ignored and one can safely use the WKB wave function 
approximately |6J| 

*i, 2 (r)~7 : ^ -IT! >< ex P ( [ [V ext (r)-^ 2 dA (61) 

{2m[V ext (r) - p]} 1/4 V V h Jr J 

with B a properly selected coefficient. The direct integration of K in eq.(^) using the above WKB wave functions is 
quite a difficult task, as already noticed in ref. || [T(| . In the simplest case the local mode of each well was assigned 
in ref. ^| the harmonic oscillator single particle ground state wave function 

and the tunneling frequency was obtained through simple integration of these Gaussian wave functions 

Sl=\n= \ [ dv$* x {r) \v ext (r) - V {2 \v - r 1)2 )l <f 2 (r) = ^-u; oe -o/«L ( 63 ) 

This result, however, is independent of the characteristic parameters of the trapped atoms, such as the number of the 
atoms, the chemical potential and the interatomic coupling constant, etc. This inadequacy is expected to be cured 
by means of the periodic instanton method presented in the following. 

Alternatively in some references the authors (29) used the external double- well potential in the form 

V ext (r) = V H (r) + V B (r) (64) 

which is created by superimposing a harmonic trap 

V H (p,z) = muj 2 (\ 2 p 2 + z 2 )/2 (65) 

and a Gaussian barrier along the z axis 

V B (z) = U exp(-^). (66) 

In their calculations, they considered a range of values for the trap frequencies u z and u) p , the barrier height Uq and 
width cr, and the condensate population N c . In the spherical or 1-dimcnsional case, one can also consider the time- 
dependent behavior by means of the variational ansatz[fl8|. Our observation is that in the vicinity of the tunneling 
region the above potential in the 1-dimensional case 



V ext (x) = -muj 2 x 2 + U exp -7— j (67) 



^muoX° + (7 exp I -—^ 2 

can be well approximated by the double-well (|68| ) considered in this paper for a variety of the geometric parameters 
of the trap. 

Now we turn to the field theory description of the GPE, i.e. the periodic instanton method which can not only 
give the correct exponential contribution of the Euclidean action but also the prefactor. The equal population case 
with N\ = N2 = N is assumed in the following calculations, due to the negligible effect of the slow time evolution 
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of the number of particles Nj (t) |p8[ . To this end we consider a scalar field problem in a 1-dimensional time plus 
1-dimensional space 



VTLLOq 2 22 

V ext (x) = -5-5- (a - x ) . 



(68) 



After a Wick's rotation t = —it the Euclidean-Lagrangian equation of motion for a finite chemical potential takes 
the form 



(69) 



The reason why we can handle a nonlinear problem by means of a linear equation of motion is that we discuss 
the tunneling behavior in the barrier region where the nonlinear interaction is negligibly small. However, there are 
obvious differences between the BEC tunneling system and the usual one-body problem, i.e. the nonlinear interaction 
contributes a finite chemical potential, which is just the integration constant on the right hand side of eq.(|6^). The 
classical turning points on both sides of the barrier can be determined by the relation V{x\ > 2) = /i as suggested in 
ref . ^] . For a noninteracting system the chemical potential approaches the ground state energy corresponding to the 
vacuum instanton case in |36|, |37]] . 

Solving this Euclidean time classical equation in the usual way J37| one obtains the periodic instanton solution in 
terms of the Jacobian elliptic function 



x c = 2xokb(k) I uiq sn (6(fc)r, k) 
where sn is a Jacobian elliptic function with modulus k, and the parameters are defined as 



(70) 



b(k) 



OA) 

2 



1 + fc 2 ' 



1 



1 



(71) 



with the imaginary time period T = 2K.{k)/b(k). The action for the half period can be calculated along the above 
instanton trajectory (70) 



where 



S, : / ' dr (\m{dx c /dT) 2 + V ext {x c ) j - IT //T/2 

T/2 



W = ~ 211(1 + u) 1 ' 2 (£(k) - uJC(ife)) . 
3 ujq 



(72) 



(73) 



Here £{k) denotes the complete elliptic integral of the second kind with modulus k. The frequency of tunneling between 
the two condensates is then given by the energy level splitting of the two lowest states, i.e. 51 = AE/fr — 2K/h = 27Z/fo 
and can be calculated by means of the path integral method as[B7| 



n = ±Ae- w ' h = ^ 
h 



2K(k' 



W 



(74) 



We emphasize here that an explicit prefactor A is included in this formula, which has been proven to be valid for the 
entire region when the chemical potential is below the barrier height. A remarkable feature of this periodic instanton 
result for the tunneling frequency is that it depends on the chemical potential fj,, or equivaientiy the number of the 
trapped atoms N through eqs. ( |49| ) or ([34]). This result will affect the conclusion about the observation of the MQST, 
as will be shown in next section. When the chemical potential approaches the top of the barrier, i.e. 



V =n 



(75) 



the periodic instanton solution ( [70| ) becomes the trivial configuration x c = with a name "sphaleron" which means 
"ready to fall" |j38| , where a type of quantum-classical transition may occur|2|]. From the eqs.(|49|) and (|60|), in the 
TFA this means 



x Q = 2R = 2a h . 



° V 2aho 



1/5 



(76) 
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where Nt = N\ + iVjj is the total number of atoms in both wells together. Therefore for a specific type of trapped 
atoms and a given double-well potential with separation x$ (the number of atomsiV^) there exists a critical number 
of atoms -ZV c i (critical separation x c \) determined by the above equation, below (above) which the tunneling process 
will give the main contribution to the tunneling amplitude. However, above this critical number of atoms or below 
this critical separation value another process, i.e. the over-barrier activation will dominate (which is definitely not 
"thermal activation" as in spin tunneling since the temperature is zero) (cf. Fig. 3). Between these two processes 
there exists a crossover. A more explicit condition for this critical number of atoms (separation between the two 



minima) can be derived beyond the TFA from eqs.(54) and (pfj), i.e 



As an example, we consider two weakly-linked condensates of Nt — 10 4 sodium atoms, confined in two symmetric 
spherical traps with frequency ujq = 100Hz as in ref. M, The critical value for x c i is x c \ = 24.58/im or more 
accurately x c \ = 25.29/zm. We note that here the height of the potential barrier is Vb = 2.21nK and the ground 
state is located at Hujq/2 = 0.38nK so that there are several energy levels beneath the barrier height. This means 
the interaction between the atoms contributes to the chemical potential, which effectively raises the classical turning 
points to a remarkably high level. Although the atoms remain in the ground state, the interaction energy is so strong 
that the vacuum instanton method can no longer be applied. We have to resort to the periodic instanton method, as 
will be shown below. 



Fig. 3 
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FIG. 3: 3 different processes through which the condensates can interchange atoms with the corresponding critical parameters 
of the boundaries. 

We now consider the "low-energy" or "non-interacting" limit, fi — > 0. As in the case of a uniform Bose gas, the 
number of atoms in the ground state can be macroscopic, i.e., of the order of the total number in one potential well, 
when the chemical potential becomes equal to the energy of the lowest state, which, in our 1 -dimensional case here, 
is fi — ► /i c = ^hujQ. The lower boundary for the chemical potential in fact implies that 



R r 



m uJo 



dho, 



(78) 



i.e. the radius of the condensate should never be less than the harmonic oscillator length <z/j . We thus have a result 
similar to that in the vacuum instanton casej37[| and the "low energy" limit here is only meaningful in this sense. 
Expanding eq. (fri|) far below the barrier height, i.e., around the modulus k — > 1, or equivalently evaluating the 
tunneling amplitude in the vacuum instanton method p7|, we obtain for the tunneling frequency 



n 



-u exp(-S c /h) 



(79) 



15 



with the Euclidean action 



S c 2 8V 2 xl 



h 3 hujo 3 af io 



(80) 



This result can be compared with that of ref. [g, eq. (p3|), which, however, gives not only a smaller exponential 
contribution 



k huj a\ o '' 



(81) 



(there is a 2/3 factor) but also an inaccurate prefactor 

A = LJ S c /h. (82) 

The source of this inaccuracy is the adoption of the too simple harmonic oscillator wave function of a single particle, 
which obviously oversimplifies the Bose-Einstein condensation tunneling problem. At least one should use the WKB 
wave function ( |6l|) in the tunneling region, and it can be shown that this corresponds to the vacuum instanton result 
we present here. For the agreement between WKB and vacuum instanton methods we refer to ref. |39| ]. 

We clarify here some features about the two-mode approximation. It is only when there are many energy levels in 
each well and the barrier height is very large that the tunneling behavior can be well defined. As already pointed out 
in ref. (§), the potential should be such that the two lowest states are closely spaced and well separated from higher 
levels of the potential, and that many-particle interactions do not significantly change this situation. This assumption 
permits a two- mode approximation to the many-body description of the system and requires that xq 3> a,ho- In the 
example considered in ref.Q the harmonic length is estimated to be au = 1.66/im for sodium atoms in a double-well 
trap with ujq — IkHz. The two-mode approximation requires that there are at least two levels beneath the barrier, 
i.e. 

1 3 

V = ^rnujlxl > -hu> , (83) 

which means the minimum separation xq should be larger than 5.76/im. For this separation eq. ( |63| ) gives the 
tunneling frequency ft ~ 0.0737-f/z, whereas the vacuum instanton result from (|79| ) is f2 ~ 2.62Hz; however, it is 
expected to reach 37% of luq in ref. M, that is 370f/z(from eq.(|63|) we know this frequency can be achieved only at a 
quite small separation xq ~ 1.6/^m). This is obviously in contradiction. To remedy this, we notice that the periodic 
instanton result eq.([74|) can reach a maximum frequency fl ~ 150Hz, which is a much higher one but for a quite 
small number of trapped atoms N c i ~ 535(Fig. 4). This can be easily understood as follows. The chemical potential 
for 5 x 10 5 atoms in the realistic experiment of the MIT group in ref. GJ, /i = 158.35nK, lies much higher than 
the barrier height Vq = 11.46nK if the two condensates are separated as close as 5.76/zm, leading to an unphysical 
parameter in eq. (|7l|) u = \J~-^ which is larger than 1. In this case the periodic instanton method cannot apply and 

the over barrier activation process will dominate over the Josephson tunneling. In the language of Josephson junction 
the two condensates seem to be connected directly with each other and no barrier exists between them. From another 
viewpoint, for the separation smaller than 5.76/im we cannot put 5 x 10 5 atoms in the double-well, i.e. the barrier 
seems no longer to exist for so many atoms and one cannot distinguish to which potential well the atoms belong and 
the oscillation is meaningless. The upper boundary for the trapped number of atoms can be estimated just to be 
N cl ~ 535 from eq.@. 

In Fig. 4 we show the TV— dependence of the tunneling frequency £1 both for a chemical potential in TFA ( |49| ) 
and beyond it (|54|). The periodic instanton results obviously lead to a rapidly growing behavior for the tunneling 



frequency when the chemical potential, i.e. the number of atoms is increased. The results of ref. |8J eq.(63) and that 
from the vacuum instanton method eq. (|79|) are also shown as horizontal lines. 



VI. OBSERVATION OF MACROSCOPIC QUANTUM SELF TRAPPING 

We are now in a position to discuss the optimal conditions for observation of the Josephson oscillation and the 
MQST effects. According to ref. for a fixed value of the initial population imbalance z(0) and phase difference 
0(0), if the parameter A exceeds a critical value A c in eq.(0), the population becomes macroscopically self-trapped 
with a nonzero average population difference (z). There are different ways in which this state can be achieved, and 
all of them correspond to the so-called MQST condition cq.([l6|). Similar results were obtained in ref.(|] for the case 
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Fig. 4 
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FIG. 4: N— dependence of the tunneling frequency Q for uj — 1kHz, xo — 5.76/im: Solid line is the result of periodic instanton 
method beyond the TFA, while dashed line that in TFA. The two horizontal lines are the results on noninteracting limit, one 
corresponds to that of ref.[^], another is from the vacuum instanton method. Inset: to = 100Hz, xo = 30fim. In both cases the 
constant results are negligibly small comparing to the N— dependent frequencies. 



that all atoms are initially localized in one well and it is concluded that the self-trapping will occur when the total 
number of the trapped atoms Nt exceeds a critical value N& determined by 

Nc2 = V ef f = — (84) 

90 Ul,2 

where V~f\ = J dr<&\ 2 (r) is the effective mode volume of each well. For this special case z(0) = 1, the critical value 
A c is shown to be A c = 2 for any initial phase difference cj>(0) and eq.(84) is obviously identical with the MQST 
condition UNt > 4K from eq.(fl6|). 

The parameters UNt and K, however, are considered as TV— independent constants in refs.[|[ Considering the 
fact that they depend actually on the number of the trapped atoms N as in our calculation above, we can refine the 
conclusions of refs. ||, g|. To access the region of self-trapping, that is, A = > A c , it is better to lower the value 
of K by making a higher barrier height Vq = ^towqITq through increasing the separation xq or the oscillation frequency 
wo 7 than to increase the number of atoms as suggested in ref. M. In fact, the quantity UNt here is proportional to 



Htf ~ N 2 / 5 as shown in eqs.(57) and (49), which means that increasing the number of atoms will not increase the 
interaction energy significantly, and at the same time the tunneling amplitude will be increased more drastically(Fig. 
4) . Thus, contrary to the observations in ref. ^ ^| , we find that the MQST will occur when the number of atoms is 
smaller (instead of larger) than a critical value AT c2 , i.e. we should decrease the number of atoms instead of increa sing 



it (Fig. 3). Inserting the values of the interaction energy eq.(p7|) or eq.(56), the chemical potential eq.([19|) or eq.(|54|) 
and the tunneling amplitude (|^) into eq.(]l6|) we may obtain this critical number of atoms for a given potential 
geometry in TFA 

8fi T F = 7A c hQ T F, (85) 

or beyond it 

8fi T F - 28C = 7A c hn b (86) 

where Qtf and Of, are the tunneling frequencies in eq.(|74|) with the chemical potential \i taking the values in TFA 
and beyond it, respectively. 

The parameters which can be adjusted are the number of atoms Nt, the oscillation frequency ojq, and the separation 
distance between the two condensates xq- Fig 3 shows the three different regions for different numbers of atoms and 
distances between the potential wells. When xq (Nt) is smaller (larger) than the critical value X2 C (N C 2), the atoms 
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Fig. 5 
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FIG. 5: Critical line for MQST effect. Solid line: results beyond TFA where the parameters take the values UNt = 8/7/itf — 4C 
and /i = fiTF + 3C/2. Dashed line: TFA results where UNt = 8/7/itf and fi = [Itf are used in the numerical simulation. 



will oscillate between these two potential wells. Once we increase the separation above (or decrease the number of 
atoms below) this critical value, the MQST will occur, i.e., most of the atoms will tend to remain in their appropriate 
wells, leading to only a small oscillation around a fixed population difference. 

We take the initial condition for the population difference to be 2(0) = 0.4 and the zero-phase case 0(0) = as an 
example. Other cases with, for example, a non-zero phase difference give rise to only a different critical parameter 
A c . For sodium atoms confined in the double-well potential with loq — 100Hz, we show numerically in Fig. 5 the 
critical line between the three different regions in TFA and beyond it. The upper region marks the self-trapping 
region, the lower the over-barrier activation. Quantum tunneling occurs only for a small range of the parameter. In 
the experiment H, the barrier was generated by an off-resonance (blue detuned) laser beam. To make our results 
more applicable to experiment we denote on the right vertical axis the corresponding barrier height in units of nK. 
We also find that the tunneling will be suppressed when the separation or the number of atoms satisfies Xo > 28/im 
or Nt > 12500 (in TFA xq > 26/xm or Nt > 12000). The crossover will occur directly between the self-trapping and 
the over-barrier regions, quite similar to the first-order transition in spin tunneling jpjjj. The self trapping region is 
very easy to access considering the easier decrease of the tunneling amplitude. 

Secondly, from the solution for z(t) in eqs.(p3|) and ( p6| ) we can define maximum amplitudes of oscillation. For 
equal populations initially, i.e. z(0) — 0, we have Hq = — cos 0(0), and the amplitude 



2 ( yjh 2 + 1 + 2A cos 0(0) - A cos 0(0) - l) ^ . (87) 



For A > 2, there is only 1 maximum, i.e., 0(0) = 7r. For A < 2, however, 0(0) = n is a minimum, while two maxima 
correspond to cos 0(0) = —A/2. Consequently from eq. (pTj) we have for A > 2 



z max = 2VA - 1/A 
and for A < 2 

W = 1 (89) 

We claim here that the result for the maximum amplitude eq.(26) of ref. |29[] is not relevant, in which the author 
erroneously took 0(0) = tt/2 as a maximum. If instead the initial population difference is nonzero, in ref. |29[ the 
author obtained the same result for the maximum amplitude as eq.([S8|) by setting 0(0) = in eq. (jl7|) . However, 
we show here this is not true for 0(0) = ir because for the ir— phase mode it is when 2(0) > 2\/A — 1/A that the 
Josephson oscillation will occur ]l0|. In this case eq.(|88|) gives the minimum amplitude instead of the maximum one. 
As already shown in Section III, the oscillation amplitude in the 0— phase mode is just the initial value of population 
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Fig. 6 




o.o ^ — 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 



z(0) 

FIG. 6: The oscillation amplitude for 0— and tt— phase modes. The slope line is the result of 0— phase mode, whereas others 
are results for different values of A, as indicated in the figure. The arrows indicate the critical value z a , which increases with 
A and finally approaches 1 for very large A. 

imbalance z(0), while in the more interesting 7r— phase mode, we plot the amplitude eq.(|35|) as a function of z(0) in 
Fig. 6 for different values of A. It is clear that the amplitudes in the n— phase mode are always larger than those in 
the 0— phase mode. The straight slope line is th e result of the 0— phase mode. For small values of A, the amplitudes 
join the slope line when z(0) > z s = y/1 — 1/A 2 , with this critical value z s increasing with A, and finally approaching 
1 for very large A. At a special value of A = 2, it is seen that the oscillation amplitude can reach 1 even for a zero 
initial population imbalance. For very large A the amplitudes in these two cases agree tangentially with each other. 
However, from Fig. 2, the condensates will be self trapped for any initial population imbalance if A > 2, making the 
Josephson tunneling unobservable. 

VII. CONCLUSIONS 

We have shown that the periodic instanton method can be used to investigate the tunneling problem in BEC 
systems at zero temperature. In particular, some deficiencies of the earlier treatments are removed. First of all, the 
tunneling amplitude and the nonlinear interaction energy between the atoms which have been taken as constants in 
refs.|^, g|, are calculated analytically in the Thomas-Fermi approximation and beyond it. The interesting features of 
the MQST effect are discussed is more detail and we observe the TV-dependence of the tunneling amplitude K and the 
self interaction energy UNt- Secondly, the N— dependence of the tunneling frequency f2 manifests a rapidly growing 
behavior when the chemical potential, i.e., the number of atoms in the trapped condensate is increased. In this sense, 
the result of ref. Q may be considered as that of a non-interacting approximation. Finnally, for the observations of 
Josephson oscillation and the self trapping effect, we suggest that one should use a small number of trapped atoms 
and change the barrier height through altering the separation xq or the oscillation frequency ujq as is evident from 
the calculation of this paper. Furthermore the ir— phase mode, that is, when the relative phase of the wave functions 
between the two condensates are opposite to each other, favors the observation of MQST, since a somewhat large A 
will bring the system out of the region of Josephson tunneling. To observe the Josephson oscillation it is better one 
adjusts instead the system into the 0— phase mode. 

Recently direct observation of an oscillating atomic current has been reported in a one-dimensional array of Joseph- 
son junctions realized with 87 Rb atomic condensate E3|. The authors verified that the BEC's dynamics on a lattice 
is governed by a discrete, nonlinear Schrodinger equation pl[ which is common to a large class of nonlinear systems. 
They used a simple variational estimate, assuming a gaussian profile for the condensate in each trap, giving the value 
of tunneling amplitude K, and a chemical potential /i ~ O.O6V0 that is much lower than the potential barrier Vq. 
This confirms our prediction that one should use a small number of trapped atoms for the observations of Josephson 
oscillation. They also observed that the wave functions, as well as K , depend on the barrier height, however, leaving 
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the analytical result unsolved. The instanton result can be extended to the investigation of the Bose-gas in such a 
periodic optical lattice fil|| , the two-component spinor condensates, or even the metastability in the case of attractive 
interaction. Work along this direction is in progress. 
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